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ABSTRACT 

We present a determination of the Hubble constant, Ho, from measurements of the Sunyaev-Zeldovich 
Effect (SZE) in an orientation-unbaised sample of seven z < 0.1 galaxy clusters. With improved X- 
ray models and a more accurate 32 GHz calibration, we obtain Ho = 64^^ ± 14 sys kms -1 Mpc" 1 for 
a standard CDM cosmology, or Ho = 66T^ ± 15 SJ(S kms -1 Mpc" 1 for a flat ACDM cosmology. In 
combination with X-ray cluster measurements and the BBN value for Ub, we find Q,m = 0.32 ± 0.05. 

Subject headings: cosmology: observations — distance scale — cosmic microwave background — dark 
matter — galaxies: clusters: individual (Coma, A399, A401, A478, A1651, A2142, 
A2256) 



1. INTRODUCTION 

The Sunyaev-Zeldovich effect (SZE) is a spectral dis- 
tortion in the Cosmic Microwave Background (CMB) due 
to inverse Compton scattering of the CMB photons off 
hot electrons; at radio frequencies this distortion is mani- 
fested as a fractional decrement in intensity of order 10~ 4 . 
For over two decades it has been known (Silk & White 
1978; Cavaliere et al. 1979) that the combination of X-ray 
and SZE observations of rich galaxy clusters, under the 
assumption of spherical symmetry, yields a direct mea- 
surement of the cosmic distance scale. Only in the past 
few years have reliable and accurate applications of this 
method become possible (e.g. Birkinshaw et al. 1991; Her- 
big et al. 1995; Carlstrom et al. 1996; Grainge et al. 1996; 
Holzapfel et al. 1997). Due to the assumption of spheri- 
cal symmetry, selection biases have been a great concern. 
To address this Myers et al. (1997) defined an X-ray flux- 
limited sample of 11 z < 0.1 clusters and began a cam- 
paign to measure distances to these clusters. This yields an 
orientation-unbiased sample, and departures from spheri- 
cal symmetry in individual clusters will average out in the 
determination of Ho from the sample as a whole. Myers 
et al. obtain a Hubble Constant of 54 ± 14 kms -1 Mpc -1 
from observations of four clusters. The accuracy of this 
result is limited by a 7% radio calibration uncertainty and 
estimated 15-30% X-ray model uncertainties for each clus- 
ter. 

In this Letter we present an improved measurement of 
Ho based on the Herbig et al. and Myers et al. re- 
sults on Coma, A478, A2142, and A2256, plus observa- 
tions of three new clusters (A399, A401, A1651) from 
their complete sample. The results we present incorpo- 
rate more accurate X-ray models (Mason & Myers 2000, 
hereafter, MM2000) and a more accurate radio calibration 
(Mason et al. 1999). We also calculate the effect of in- 
trinsic CMB anisotropies on our result and find them to 



be a limiting factor for a sample of this size. Throughout 
this paper we use Ho = 100/ikms -1 Mpc -1 , and consider 
two cosmologies: SCDM (tt M = l,fi A = 0), and ACDM 
(f2j\/ = 0.3, ft\ — 0.7). Unless otherwise stated, we assume 
the SCDM cosmology. 

2. THE SUNYAEV-ZELDOVICH EFFECT, CLUSTER 
SAMPLE, AND X-RAY MODELS 

2.1. The Sunyaev-Zeldovich Effect 

The observed fractional change in antenna temperature 
induced by a cloud of electrons with a beam-averaged 
Compton y— parameter y f, s is (Sunyaev & Zeldovich 1980) 

= (11^2 - 4] Xrel X Vcba- (1) 

cmb I, J- ) 

Here x = hv/kTcmb, and Xrel 1S a relativistic correction 
factor which we compute using the analytic expression of 
Sazonov & Sunyaev (1998). We assume zero peculiar ve- 
locity, which introduces a ~ 5% uncertainty but will aver- 
age out over the sample. Along a given line of sight, the 
Compton y— parameter is 

V=— 2 r (2) 

m e cr 

where r is the inverse Compton optical depth. At 
32.0 GHz (x = 0.563), Eq. 1 takes the form 

= -1.897 xy obsX ^. (3) 

crab 

We follow the procedure of Myers et al. in correcting 
the observed Compton y— parameters for relativistic ef- 
fects rather than the models. 

Given a model for the cluster density and temperature 
profiles from X-ray observations, we predict the SZE decre- 
ment in terms of the Compton y— parameter averaged over 
the telescope beam and switching patterns: 

Upred^-^- / dVLy mod el(P) Rn,sw{®)- (4) 
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Here flbeam is the solid angle of the telescope beam, y m odei 
is the predicted y along a given line of sight, and Rn,sw(^) 
is the normalized beam response, including the effects of 
switching. It is easily shown that y pre d, as determined 
from X-ray observations and Eq. 4, is proportional to 
hT 1 !" 1 , so that if y a b s = q and y pre d — phT x l 2 , the Hubble 
constant is given by 



Note that the t sw values presented in MM2000 are simply 
the model inverse Compton optical depths convolved with 
the telescope beam switching pattern as per Eq. 4. For an 
isothermal cluster, y prec i = t sw kT e /m e c 2 . 

2.2. Cluster Sample 

In MM2000 we define a larger, X-ray flux-limited clus- 
ter sample which expands upon that presented by My- 
ers et al. . This is a 90% volume-complete sample 
selected from the XBAC catalog (Ebeling et al. 1996) 
with z < 0.1 and F x > 1.0 x 10~ n ergcm~ 2 s _1 . At 
z = 0.1, the volume-completeness criterion corresponds to 
L x > 1.13 x 10 44 ft- 2 ergs- 1 (0.1 - 2.4 kcV). The resulting 
set of 31 clusters contains the Myers et al. sample. While 
the 7 clusters for which we present measurements here are 
all members of the smaller Myers et al. sample, this work 
is part of an ongoing project to survey distances to the ob- 
jects in our expanded sample. MM2000 discusses in detail 
the completeness of this XBAC-derived catalog. 

2.3. X-Ray Models 

MM2000 also presents X-ray models for the 22 clus- 
ters in our sample which had public ROSAT PSPC data 
as of May 1999. The primary focus of this analysis was 
to quantify the uncertainty in the beam-convolved inverse 
Compton optical depth, t sw , a quantity which we find is 
robustly constrained by the ROSAT data. We adopt the 
Markevitch et al. (1998) measurements of the cluster gas 
temperature T e , which — except for the Coma cluster — we 
assume to be isothermal. For the Coma cluster, following 
Hughes ct al. (1988), we adopt a hybrid model which is 
isothermal inside a radius of 500 h^ 1 kpc, outside of which 
the temperature follows the density profile with an adi- 
abatic index of 7 = 1.5. For A478 and A2142, we use 
the fits which have been corrected for the spectral bias 
induced by the cooling flow emission. Table 1 summa- 
rizes our cluster models and the resulting beam-averaged 
Compton y— values, y pre d- In § 4 we discuss the impact on 
Ho of our choices for the cluster models. 

3. SZE OBSERVATIONS 

In this section we report on observations of 7 clusters 
from our sample. Two of these, A399 and A401, are in 
close proximity. The effects of this are included in our 
model predictions (see § 4). 

3.1. New Observations 

For the observations reported here, the OVRO 5.5- 
meter telescope (5-m) was outfitted with a HEMT re- 
ceiver having a center frequency of 32 GHz and a band- 
width of 6.5 GHz. The HEMT input is Dicke switched 



every millisecond between two ambient-temperature cor- 
rugated feeds which give rise to two 7'. 35 FWHM Gaus- 
sian beams [each having a main-beam volume of Qbeam = 
(5.21 ±0.03) x 10" 6 Sr] separated by 22'. 16 in azimuth. To 
remove systematic effects due primarily to atmosphere and 
ground, we employ a triple-differencing technique. Two 
levels of differencing are provided by the Dicke switching 
between the target (ON) position and a reference (REF) 
position, and and by nodding the telescope in azimuth at a 
rate of ~ 0.1 Hz. The third level is provided by observing 
blank leading (LEAD) and trailing (TRAIL) fields over the 
same range of azimuth and elevation as the MAIN field, so 
that ground-based signals are cancelled in the difference. 
The cluster signal is given by the MAIN field signal minus 
the average of the LEAD and TRAIL signals, or MLT. 
The LEAD-TRAIL field difference (LTD) provides a diag- 
nosis of possible residual signals due, e.g., to unsubtracted 
ground emission or intrinsic CMB fluctuations. 

Flux density calibration is accomplished with reference 
to Cas A using an epoch 1998 flux of 194 ± 4 Jy (Mason 
et al. 1999), with the secular variation of Cas A modeled 
as per Baars et al. (1977). Data are edited by automatic 
filtering algorithms that use the scatter of the data and 
their standard deviations as rejection criteria. More dis- 
cussion of our observing strategy and data filters is given 
in Readhead et al. (1989) and Myers et al. (1997); more 
details on the instrument are in Leitch et al. (2000). 

3.1.1. Discrete Source Removal 

The brightest contaminating discrete sources in A399 
and A1651, selected from 87GB (Gregory & Condon 1991), 
were monitored at 18.5 GHz with the OVRO 40-meter tele- 
scope in November and December of 1996, concurrent with 
cluster SZE observations. In the Fall of 1997 we used the 
OVRO 40-m telescope to survey all NVSS (Condon et al. 
1998) sources in these fields and those for A401 which, 
assuming a flat spectrum, would give a peak signal at 32 
GHz > 16 [iK. We also observed all sources within 12' of 
any of the field centers and having 1.4 GHz flux densities 
> 50 mjy. Sources which showed indications of variability 
were extrapolated to 32 GHz assuming a flat spectrum; 
otherwise the two-point spectral index between 1.4 and 32 
GHz was used. The overall corrections in are given 
for each cluster in § 3.1.2. 

3.1.2. A399, A401, and A1651 

Observations of A399 were taken from October 1996 
through March 1997. For these observations the LEAD 
and TRAIL fields were separated from the MAIN by 26 
minutes of Right Ascenscion. After statistical filters and 
weighting of the data, we acquired 40 hours of total in- 
tegration time including LEAD and TRAIL fields, and 
time spent on the REF beams. To account for source 
contamination 24 ± 5/uK was subtracted from the ob- 
served signal. We then determine a decrement of AT & S = 
-164 ± 21 (MLT). The LEAD-TRAIL difference is 
ATltd = 15 ± 21 /j,K. The observed AT corresponds to 
2/ofcs = (3.28 ± 0.42) x 10~ 5 , including the relativistic cor- 
rection Xrel — 1-027. 

A1651 was also observed during this period, with LEAD 
and TRAIL field separations of 20 min 30 sec from the MAIN 
field; the effective total integration time was 25 hours. We 
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Table 1 
Cluster Model Parameters 



Cluster 


z 


n 


a 
P 




kT e 








(') 




(lO" 3 h 1 / 2 cm- 3 ) (keV) 


(io- 3 h- 1 ' 2 ) 


A399 


0.0715 


4.33 ± 0.45 


0.742 ± 0.042 


3.23 ±0.18 


7.0 ±0.2 


2.50 ±0.13 


A401 


0.0748 


2.26 ±0.41 


0.636 ± 0.047 


7.90 ±0.81 


8.0 ±0.2 


3T7±0.18 


A478 


0.0900 


TOO ±0.15 


0.638 ±0.014 


27.81 ±9.7 


8.4 ±0.7 


3.68 ±0.15 


A1651 


0.0825 


2.16 ±0.36 


0.712 ±0.036 


7.14 ±3.20 


6.1 ±0.2 


2.44 ±0.11 


Coma 


0.0232 


9.32 ±0.10 


0.670 ±0.003 


4.52 ±0.04 


9.1 ±0.4 


2.76 ±0.16 


A2142 


0.0899 


1.60 ±0.12 


0.635 ±0.012 


14.95 ± TO 


9.7 ±0.8 


4.28 ±0.18 


A2256 


0.0601 


5.49 ±0.21 


0.847 ±0.024 


4.08 ± 0.08 


6.6 ±0.2 


3.19 ±0.16 


NOTE: 


Tsw, do, (3, and n eo are 


from MM2000; T e values are 


from Markevitch et al. (1998); 



redshifts are from Struble & Rood (1991). Uncertainties arc la. 



subtract 5 ± 1 /x.K from the observed decrement to cor- 
rect for discrete sources. These are mostly in the REF 
beam of the LEAD field. With the source correction we 
determine a decrement AT obs = -247±30 /xK (MLT). The 
LTD is AT LTD = -92 ± 36 /xK. The Compton-jy parame- 
ter is y bs — (4.88 ±0.59) x 10~ 5 including the relativistic 
correction \rei = 1.023. 

A401 was observed from October 1997 to March 1998, 
giving an effective total integration time of 43 hours. The 
LEAD and TRAIL field separations were 16 min 36 scc . Dis- 
crete sources contribute 16 ± 11 /iK. The source-corrected 
decrement is AT obs = -338 ± 20 /iK (MLT); the LTD is 
AT LTD = 133 ± 24 /iK. This gives (with X rei = 1-031) 
y obs = (6.33 ±0.43) x TO" 5 . 

The statistically significant LTD's for A401 and A1651 
are discussed in § 4. 

3.2. Coma Cluster 

Herbig et al. (1995) used the OVRO 5-m telescope to 
determine the SZE decrement towards the Coma cluster 
(Abell 1656) during the observing seasons of 1992 and 
1993, and found a Rayleigh- Jeans temperature decrement 
of AT bs = —302 ± 48itK. These data were calibrated 
relative to DR21 assuming Sdr21 = 18.24 ± 0.55 Jy, and 
referenced to the telescope main beam assuming flseam = 
(5.16 ± 0.15) x 10~ 6 Sr. From the Mason et al. flux den- 
sity scale and Xrei — 1.035, we determine a beam-averaged 
Compton y— parameter of y sw = (6.38 ± 1.01) x 10~ 5 . 

3.3. A478, A2142, and A2256 

A478, A2142, and A2256 were observed by Myers et al. 
(1997) from July 1993 to March 1994 resulting in decre- 
ments of -375 ± 28, -437 ± 25, and -243 ± 29 /xK re- 
spectively; these temperature are referred to the main 
beam, and have been corrected for point-source contam- 
ination. These data were calibrated using a brightness 
temperature of Tj = 144 ±8 K for Jupiter. After in- 
cluding the correction for the main beam area [Myers et 
al. use (5.12 ± 0.14) x 10~ 6 Sr], we find an overall correc- 
tion of / = 1.037 for the Myers et al. data. The stated 
decrements yield switched, beam-averaged Compton-j/ pa- 
rameters of (7.25 ± 0.54) x 10" 5 , (8.44 ± 0.48) x 10~ 5 , 
and (4.70 ± 0.56) x 10~ 5 , respectively. We apply rela- 
tivistic corrections (xrei = 1-033 , 1.037 , and 1.026) to- 
gether with our calibration correction, and find Compton 



y-values of (7.77 ± 0.58) x 10" 5 , (9.10 ± 0.52) x 10" 5 , 
and (5.00 ±0.60) x 10~ 5 . 

4. INTERPRETATION 

For SZE observations at our frequencies and angular 
scales intrinsic CMB anisotropies are a significant source of 
uncertainty. We use the RING5M measurements (Leitch 
et al. 2000) to determine the impact of this on our H 
results. In measurements taken with the OVRO 5-m tele- 
scope at 32 GHz, Leitch et al. determined an RMS tem- 
perature difference on 22' scales of 5T 2 2> — 79 J ^K- In 
order to determine the degree to which parallactic angle 
averaging during a track on a cluster reduces this, we 
generated 10 3 realizations of 4deg 2 patches of sky from 
a ACDM power spectrum. Each realization was convolved 
with the 5-m main beam and the switching pattern char- 
acteristic of a typical SZE observation. We find (including 
the LEAD/TRAIL differencing) a residual CMB signal of 
81.6% x 6T 22 >, or 64 /iK (a y = 1.24 x 10" 5 ) assuming the 
RING5M power level. The RMS of LTD predicted by this 
analysis is 90 /iK, close to the observed RMS of 79 /iK in 
our 7 clusters. On this basis the non-zero LTD's seen in 
A1651 and A401 arc not unexpected. 

While this noise level is not small compared to the sig- 
nal in fainter clusters such as A399, we leave these clusters 
in our sample to avoid an orientation bias and perform a 
Maximum Likelihood analysis to determine H . Assume 
that we have N observed Compton- y paramaters q i: and 
N predicted y-parameters (as per Eq. 5) pi determined 
from X-ray observations. The SZE observations are cor- 
rupted with random Gaussian noise o q ^ (including ther- 
mal noise and point-source subtraction uncertainties) and 
the predictions are corrupted with random Gaussian noise 
cr Pi i; the RMS residual CMB signal is a cm \, = 64 /xK. In 
addition there are N "true" predictions pi which are the 
predictions which we would have in the absence of errors 
in the X-ray models. Then we can determine h by mini- 
mizing 

2 _ ^ fa - fe~ 1/2 f t ) 2 (Pi -Pi) 2 
i=\ q-i + crab "p.i 

where x 2 = —2lnL. The unconstrained pi must also be 
minimized in this process, but these can be projected out 
analytically. Setting the derivatives of x 2 with respect to 
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the pi equal to zero, we find 

~ _ qi + fih 1/2 pj 

with fi — (a 2 ; t + cr^ mh )/CTp ^. This leaves a 1-dimensional 
parameter space. The value of h which minimizes Eq. 6 
is the most probable value; two-sided 68% confidence in- 
tervals are determined about this value by integrating the 
likelihood function L. Wc impose the prior 0.01 < h < 5; 
our upper bound on H for A399 is slightly affected by 
our choice of prior, but other results are not. Monte Carlo 
tests with the estimator of Eq. 6 show that within the 
range of power levels allowed by RING5M, our result is 
not affected by the uncertainty in <7 cm b- 

Table 2 shows the measured and predicted values of the 
Compton y— parameter for our 7 clusters along with the 
Maximum-Likelihood values for and 68% confidence lim- 
its on Ho. A399 and A401 are in close spatial proximity 
with the result that the observed decrement on either is 
reduced somewhat due to signal from the other in the ref- 
erence beam; this is accounted for in the stated values of 
Upred- Averaged over a track, we find a decrease in the 
predicted y for A399 of (0.15 ± 0.05) x 10 -5 /i -1 / 2 and for 
A401 of (0.13 ± 0.05) x 10 -5 /i -1 / 2 . For an SCDM cos- 
mology, the sample average is 64 kms -1 Mpc -1 . The 
scatter of the H values gives an error in the mean of 
only 8 kms -1 Mpc -1 , which is less than that given by 
our Maximum Likelihood analysis. The minimum of \ 2 
corresponds to xl — 0-32 for v — 6; there is a ~ 5% 
chance of obtaining a xt value this low by chance for 
six degrees of freedom. We therefore think it likely that 
the scatter in these 7 clusters is fortuitously small. Since 
our observing and modelling errors have been well quan- 
tified, the Maximum Likelihood method gives a more re- 
liable estimate of the uncertainty than the data scatter 
for a small sample such as ours. Our result is not sensi- 
tive to the model we choose for the cluster temperature 
profile: adopting hybrid models for all 5 clusters reduces 
the sample average by only 4%. Although more extreme 
models for the temperature profiles would have a larger 
effect on our result, such models are not motivated by 
current X-ray analyses (Irwin, Bregman & Evrard 1999; 
White 2000). For a ACDM cosmology, the sample average 
is 66^11 km s _1 Mpc -1 . The calibration uncertainties are 
3% (radio) and 8% (X-ray), and we estimate a 10% uncer- 
tainty in the SZE predictions due to the possibility of sub- 
structure and non-isothermality in the ICM. Altogether we 
have a systematic error budget of ±14 sys kms -1 Mpc -1 , or 
±15 sas kms -1 Mpc -1 for ACDM. 

The increase in H relative to the Myers et al. result 
is due primarily to differences in the updated X-ray mod- 
els. As a check on our ROSAT-derived density models, we 
have compared the baryonic masses that we obtain inside 
500 h- kpc to those obtained by Mohr et al. (1999) for 20 



of the 22 clusters in common to our analyses: we find that 
the mean mass ratio is equal to unity at better than 1% 
accuracy. See MM2000 for more discussion of this anal- 
ysis. Also the T e values we have used tend to be higher 
than the Myers et al. values due to the correction for the 
spectral bias of the cooling flows. Using the Myers et al. 
temperatures increases the fractional scatter in the Ho val- 
ues by over a factor of two and decreases the sample mean 
to 55 kms -1 Mpc -1 . The fact that the uncertainty in our 
result is comparable to that of Myers et al. — in spite 
of improved X-ray models, better calibration, and 3 more 
clusters — is due to our inclusion of intrinsic anisotropy in 
the analysis. 

5. CONCLUSIONS 

We have presented a determination of Ho resulting from 
measuring the SZE in 5 clusters from an unbiased, low- 
rcdshift sample. For an SCDM cosmology we find Ho = 
64±]$ ± 14 sy;s kms -1 Mpc -1 , while for a ACDM cosmol- 
ogy we obtain H = 66+^ ± 15 sys kms -1 Mpc -1 . This 
result is in good agreement with other recent measure- 
ments. Madore et al. (1999), using Cepheid variables as 
distance indicators, find 72±5±7 sys kms -1 Mpc -1 . Grav- 
itational lens time delays give results consistent with these 
(e.g. Fassnacht et al. 1999). 

Our value of Ho implies an age for the universe ranging 
from (10.2 ± 3.2) x 10 9 years for an SCDM universe to 
(14.2 ± 4.5) x 10 9 years for ACDM. These are both con- 
sistent with recent age determinations from main-sequence 
fitting, which give ages of (> 12±1) x 10 9 years (Chaboyer 
et al. 1998). In the X-ray cluster analysis of MM2000 we 
find a mean ICM mass fraction ficM = Mb a ry/M t ot = 
(7.02 ± 0.28) /i -3 / 2 x 10 -2 (within R 500 ~ l/i -1 Mpc) in 
a sample of 22 nearbly clusters. We combine this result 
with the Big Bang Nucleosynthesis constraint Sl^/i 2 = 
0.019 ± 0.001 (Buries & Tytler 1998) and our value for 
Ho to find a total matter density Am = 0.32 ± 0.05, argu- 
ing against SCDM cosmologies. 
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Table 2 
h results on 5 clusters 



Cluster y obs y pred H 

(1(T 5 ) (lO- 5 ^- 1 / 2 ) (kms^Mpc- 1 ) 
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3.27 ±0.20 


102±^ 
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48t 2 * 


A478 
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6118 
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Coma 
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62+ 49 

DZ -24 
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79±H 
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67ti 


SAMPLE 
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